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Abstract. We report new near-infrared, long-baseline interferometric observations at the AU scale of the pre- 
main-sequence star FU Orionis with the PTI, IOTA and VLTI interferometers. This young stellar object has been 
observed on 42 nights over a period of 6 years from 1998 to 2003. We have obtained 287 independent measurements 
of the fringe visibility with 6 different baselines ranging from 20 to 110 meters in length, in the H and K bands. Our 
data resolves FU Ori at the AU scale, and provides new constraints at shorter baselines and shorter wavelengths. 
Our extensive (u, nj-plane coverage, coupled with the published spectral energy distribution data, allows us to 
test the accretion disk scenario. We hud that the most probable explanation for these observations is that FU Ori 
hosts an active accretion disk whose temperature law is consistent with standard models and with an accretion 
rate of M = (6.3 ± 0.6) x lt)~^[Mi,/M q)~^ Mq/yi. We are able to constrain the geometry of the disk, including 
an inclination of 551® deg and a position angle of 47j[l deg. In addition, a 10 percent peak-to-peak oscillation is 
detected in the data (at the two-sigma level) from the longest baselines, which we interpret as a possible disk hot¬ 
spot or companion. The still somewhat limited [u, v) sampling and substantial measurement uncertainty prevent 
us from constraining the location of the spot with confidence, since many solutions yield a statistically acceptable 
ht. However, the oscillation in our best data set is best explained with an unresolved spot located at a projected 
distance of 10 ± 1 AU at the 130 ± 1 deg position angle and with a magnitude difference of AK « 3.9 ± 0.2 and 
AH ~ 3.6 ±0.2 mag moving away from the center at a rate of 1.2 ±0.6 AUyr“^. Although this bright spot on the 
surface of the disk could be tracing some thermal instabilities in the disk, we propose to interpret this spot as the 
signature of a companion of the central FU Ori system on an extremely eccentric orbit. We speculate that the 
close encounter of this putative companion and the central star could be the explanation of the initial photometric 
rise of the luminosity of this object. 

Key words. Stars: pre-main-sequence, late-type, individual: FU Ori - Planetary systems: protoplanetary disks - 
Infrared: stars - Accretion, accretion disks - technique: interferometric 


1. Introduction 

FU Orionis is a variable stellar system representative of 
a small class of pre-main-sequence (PMS) objects named 
FUor s llAmbartsumianI fl971t lAmbartsumian fc Mirzoi^ 
Il982h . These stars show a rise of their luminosity over a 
timescale of a few hundred years. Some FUors, observed 
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before their flare-ups, display the spectral characteristics 
of a typical T Tauri star (TTS). It is now widely accepted 
that most TTS go through this type of very short FUor- 
like outburst phase, possibly several times during the 
early stages of their existence ilHerbidll977t iKenvon et alJ 
l2nn(1i . 

FUors display other interesting observational peculiar¬ 
ities: an infrared excess that is far larger than the one 
found in T Tauri stars, and very broad absorption lines 
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that appear to be double. The exact nature of FUors prop¬ 
erties is still controversial, with essentially two competing 
explanations: 


1. Rapidly Rotating G Sup ergiant: IPetrov fc Herbid 
ll1992^ ; iHerbig et al.l ll200,‘lli have shown that an FUori 
optical spectrum can be reproduced by a G-type super¬ 
giant chromosphere overlaid with a r i sing, cooler ab¬ 
sorbing shell. As proposed bv Ibarsoi] |Ii 98(TI1 . the star 
is rotating near breakup, and could be responsible for 
the flare-up. 

2. Accretion- disk Instability: iHartmann fc KenvonI 
ill 98,^17^ have proposed that the flare-up is a phe¬ 
nomenon not of the star itself but rather the result of a 
major increase in the surface brightness of the circum¬ 
stellar accretion disk surrounding a young T Tau-type 
star. Since the idea was first introduced, theories of in¬ 
stabilities int rinsic to such an acc re tion disk have bee n 
examined by IClarke et al.1 lll99(il : iBell fc Lh] lll994ll : 

an d iKlev fc* T;inl l|l 99*^ . A related scenario 
is the passage of a companion star close to the primary 
star, which may cause an increase in the accret ion rate 
llBonnell fc Bastienlll992l: IClarke fc Sveilll9^ . 


IHartmann fc KenvonI ^ll99(^^ have three main argu¬ 
ments to explain the origin of spectroscopic properties 
of FUors with the disk accretion model. First, the sep¬ 
aration of the peaks in the double absorption lines de¬ 
creases at longer wavelengths. Such behavior is expected 
in a Keplerian accretion disk, since the outer parts ro¬ 
tate slower than the inner ones. Second, the presence of 
CO bands requires environmental conditions cooler than 
the stellar surface. Such conditions can be found in the 
inner regions of the disk. Third, as is the case for other 
young stellar objects (YSOs), the infrared excess produced 
by a disk can explain the shape of the observed FUor 
infrared spectral e nergy distributions (SEDs). However 
iHerbig et alJ ll200,8ll claim that “the spectroscopic proper¬ 
ties of FUors that have been argued as proof of the accre¬ 
tion disk picture [...] are found in a number of much older 
high-luminosity stars where the re is no reason to thin k 
that an accretion disk is present” . IHartmann et al.1 1 2flfl4l . 
based on additional optical spectroscopic observations of 
FUors, point out that the chromospheric emission reversal 
of lines with very different strengths and excitation poten¬ 
tials required in the framework of a rapidly rotating star, 
as well as the decrease of the rotational broadening with 
wavelength, do not favor the chromosphere model. 

High angular resolution observations of these objects 
can bring new observational constraints and help to solve 
this controversy. FU Ori is located in Orion [d = 450 pc; 
1 AU correspond s to 2.2 mas . ) Rec ent observations with 
adaptive systems FWang et al showed that FU Ori 

has a companion located 0.5arcsec south with a magni¬ 
tude difference of about 4 magnitudes in the near infrared. 
This companion has been found to also be a young star 
iReipiirth fc Aspiiil ||20 m3). 

The first high angular resultion observations of PMS 
stars with an infrared interferometer were carried out on 


FU Ori itself llMalbet et al.lll9^ , hereafter paper I), 
one of the brightest YSOs, using the Palomar Testbed 
Interferometer (PTI). Unfortunately the limited number 
of baselines was not sufficient t o definitively conclu de on 
the exact nature of FU Orionis. iBerger et al.1 l|200nll pub¬ 
lished additional visibility data on FU Orionis obtained 
with PTI and the Infrared and Optical Telescope Array 
(IOTA). 

In this paper, we present the results of a 6 year cam¬ 
paign on FU Ori with the most advanced infrared interfer¬ 
ometers around the world, including the recently available 
Very Large Telescope Interferometer (VLTI). Section [21 
summarizes the observations; Sect. O describes briefly the 
data processing and Sect. 0 presents the results. We ana¬ 
lyze and interpret these new observations in Sect. Eland 
we discuss them in the framework of disk models and star 
formation theories in Sect. El 


2. Observations 

Observations were carried out with three long-baseline 
near-infrared interferometers: IOTA, PTI, and VLTI. 

IOTA, located on Mount Hopkins, was at the time 
of the observ ations a two-telescope stellar interferometer 
llTra,iiblll998ll . whose L-shape multiple baseline configu¬ 
ration allows observations with baselines ranging from 5 
to 38 m, with limiting magnitudes of AT ~ 6,11 ~ 6. 
PTI is a t hree-telescope interfe rometer located on Mount 
Palomar llColavita et al.lli999^ with limiting magnitudes 
of AT ~ 6.5, iL ~ 5 and whose 3 baselines can be operated 
separately. The VLTI is an interferometer that at the time 
of the observations could combine two apertures among 
four 8m unit telescopes (UTs) or two 0.4m siderostats 
llGlindemann et ^l2003ll . The VLTI was used with tip- 
tilt correcti on on UTl and UT3 and the VINCI focal 
instrument llKervella et all l200fl l2003ll which is a fiber- 
filtered instrument capable of combining two beams in the 
AT band. 

FU Orionis was observed during five observing runs in 
1998 with IOTA, in 1998, 1999, 2000, 2003 with PTI, and 
in 2002 with the VLTI. We spent a total of 42 nights on 
this project distributed into 29 nights on PTI, 11 on IOTA 
and 2 on VLTI. The distribution of these observations 
with time and baselines is summarized in Tabled A total 
of 287 observations have been performed leading to 287 
fringe visibility points. 

Table 121 summarizes the main characteristics of our ob¬ 
servations: dates, baselines, filters and calibrators. 


3. Data processing 

The data processing involves two main steps: 

1. The extraction of raw visibilities from recorded data, 
a procedure that is specific to each interferometer. 

2. The data calibration, i.e. dividing raw visibilities by 
an estimate of the instrumental transfer function, a 
procedure that is common to all data sets under the 
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Table 2. Logs of FU Ori observations. 


Observing period 

Baselines 

Filters 

Calibrators 

1998-NOV-14 

1998-NOV-27 

PTI 110m (NS) 

K 

HD 42807, HD 37147, HD 38529, HD 32923 

1998-Dec-13 

1998-Dec-26 

IOTA 21m (S15N15), 38m (S15N35) 

K, H 

HD 42807, HD 37147, HD 38529, HD 31295 

1999-NOV-23 

1999-Dec-Ol 

PTI 110m (NS) 

K, H 

HD 42807, HD 37147, HD 38529, HD 32923 

2OOO-N0V-I8 

2OOO-N0V-27 

PTI 110m (NS), 85m (NW) 

K, H 

HD 42807, HD 37147, HD 38529, HD 32923 

2002-Oct-28 

2000-Oct-29 

VLTI/VINCI 100m (UT1-UT3) 

K 

HD 42807 

2OO3-N0V-I9 

2OO3-N0V-27 

PTI 110m (NS), 87m (SW) 

K 

HD 42807 

Angular diameters estimated from Hipparcos catalog in mas: 

HD 42807: 0.5 ± 0.1, HD 37147: 0.3 ± 0.2, HD 38529: 0.4 ± 0.3, HD 31295: 

0.5 ± 0.2, HD 32923: 0.75 ± 0.30. 


Table 1. Distribution of observations over the campaign. 





# nights 


# 

Baseline 

Band 

98 

99 00 02 

03 

Data 

IOTA/S15N15 

H 

3 



44 

IOTA/S15N35 

H 

5 



39 

PTI/NS 

H 


2 


13 

PTI/NW 

H 


2 


4 

total 



12 


100 

IOTA/S15N15 

K 

1 



4 

IOTA/S15N35 

K 

2 



18 

PTI/NS 

K 

11 

4 1 

2 

109 

PTI/NW 

K 


4 


32 

PTI/SW 

K 



3 

10 

VLTI/U1-U3 

K 


2 


14 

total 



30 


187 


assumption that the raw visibilities have had all in¬ 
strumental biases corrected by the previous step. 


The IOTA interferometer temporally encodes fringes. 
We chose a quadratic estimator similar to the one used 
for FLUOR llCoude Du Foresto et ahlUbQTll with no pho¬ 
tometric calibration signal. The two interferograms simul¬ 
taneously recorded at the output of the instrument are 
subtracted to build a single interferogram with reduced 
photometric contamination^. Visibilities are computed by 
estimating the energy contained at the fringe position in 
the spectral power density. We compute one average visi¬ 
bility from each batch of 500 interferograms and the stan¬ 
dard deviation provides the error estimate on this mea¬ 
surement. 

The VINCI interferograms are, much like the IOTA 
case, also produced by a temporal modulation of the 
optical path difference. The processing applie d to this 
signal is described in iKervella et al J (|20^, and is 
based on a strategy compar able to that defined by 
ICoude Du Foresto et all lll997l). with the optional use of 
wavelet transforms - as opposed to Fourier transforms - 
for fringe detection and power spectrum computation. As 
in the case of the IOTA data, we chose to use the visibil¬ 
ities computed from the classical Fourier power spectrum 
method. Unlike the IOTA case, we used the available pho¬ 
tometric channels to compensate for coupling efficiency 
variations. Each visibility measurement results from the 
averaging of a 500-interfero CTam set of observations a nd 
the bootstrapping method l|Efron fc Tibshir;milll99,3ll is 
applied to derive its associated internal error. 


3.1. Extracting raw visibilities 

PTI data processing is described bv IColavit^ lll999ll . The 
visibility is based upon the ABCD algorithm which com¬ 
putes its estimate from 4 different phase quadratures of 
one fringe. Out of the four fringe visibility estimation al¬ 
gorithms available in the standard PTI data processing 
pipeline we chose the incoherent spectral estimator; the 
spectrometer has a higher instrumental transfer function 
(due to the spatial filtering effect of the single mode fiber), 
and the incoherent estimator allows one to average visibil¬ 
ities over the entire H or K bandpass without introducing 
biases due to atmospheric piston (cf. paper I). 


3.2. Data external calibration 

After the computation of raw visibilities for each interfer¬ 
ometer, the data reduction process enters the calibration 
path common to all the interferometers. The key point is 
to estimate the system transfer visibility for each obser¬ 
vation of the source by observing calibration sources for 
which one can estimate the fringe visibility that would 
be measured with an unresolved source. These calibration 

^ Since the fringes are in phase opposition, the energy in the 
spectral density distribution is maximized at the fringe position 
and the photometric energy at lower frequencies is reduced. 
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sources were chosen to minimize time- and sky-dependent 
variations, hence were close on the sky, and observed 
within 30 minutes of the science target. The calibrators 
are corrected for their intrinsic visibility loss due to their 
apparent diameter (inferred from Hipparcos distance and 
size based on spectral type). We estimate the system 
transfer visibility at the time of the target observations 
by interpolating the calibrator diameter-corrected visibil¬ 
ities. The weight of each calibrator measurement is given 
by the inverse of the time delay between the calibrator 
measurement and the target measurement. Division of the 
observed target visibility with the interpolated system vis¬ 
ibility leads to the determination of the calibrated visibil¬ 
ity and its associated error. 

4. Results 

The observations with all interferometers covered a range 
of hour angles large enough to allow the projected base¬ 
lines on the sky to rotate over a significant angle. Figure 
n displays the resulting (u, v) plane coverage of our obser¬ 
vations with the six baselines: 21m (IOTA-S15N15), 38m 
(IOTA-S15N35), 85m (PTI-NW), 87m (PTI-SW), 100m 
(VLTI-U1U3) and 110m (PTI-NS). 

Figure El shows our results for the six baselines in the 
H and K bands (left and right columns respectively). All 
panels show calibrated squared visibilities as a function of 
hour angle. The main raw results are: 


Table 3. Averaged square visibilities for the 6 baselines. 


baseline 

A 

Bp (m) 

e (deg) 

lu^l 

IOTA/S15N15 

H 

20.2 ±0.5 

104.5 ±9.5 

1.00 ±0.13 

IOTA/S15N35 

H 

35.0 ±0.7 

95.7 ± 13.7 

0.91 ±0.14 

PTI/NS 

H 

103.7 ±2.2 

61.8 ±3.5 

0.83 ±0.04 

PTI/NW 

H 

85.6 ±1.5 

14.1 ±0.8 

0.79 ± 0.05 

IOTA/S15N15 

K 

20.8 ±0.3 

114.9 ±3.4 

0.90 ±0.13 

IOTA/S15N35 

K 

34.4 ± 0.3 

91.4 ±7.8 

0.98 ±0.09 

PTI/NS 

K 

102.7 ±3.8 

63.6 ±6.5 

0.72 ± 0.08 

PTI/NW 

K 

84.2 ± 1.8 

14.0 ± 1.8 

0.79 ± 0.05 

PTI/SW 

K 

82.5 ±2.9 

-55.6 ±4.9 

0.82 ±0.08 

VLTI/U1-U3 

K 

89.6 ± 7.2 

54.2 ±5.3 

0.87 ±0.05 


oscillation as function of the spatial frequency. This can 
be explained by: 

— the presence of a source which is marginally resolved 
and which dominates the overall flux, 

— and by the presence of a second, fainter source. The 
fainter source can produce the oscillatory feature if it 
is located at a distance larger than the typical size of 
the compact centered structure. 


— We confirm the 1997 result (paper I): FU Orionis is 
resolved at 110m in the K band with the same visibility 
range (PTI-NS data) of 0.72 ± 0.07 (cf. paper I). 

— FU Orionis is also resolved at about the same visibility 
level with the 85m PTI-NW, 87m PTI-SW, and 100m 
VLTI-U1-U3 baselines. 

— FU Orionis is unresolved at short baselines (IOTA 
data). 

— We obtain color information between 1.6/rm (H-band) 
and 2.15/im (K-band). 

— We detect a peak-to-peak oscillation of 10% at 110m in 
the K band and more marginally in the H band (this 
oscillation appears better in binned data of PTI/NS 
displayed in the right panel of Fig. El). 

5. Interpretation of FU Ori observations 

Interpretation of the observations performed with infrared 
long-baseline interferometers is not obvious since we have 
no direct information about the brightness distribution 
of the object. However, the measured visibilities allow us 
to set constraints on current models. In paper I, we had 
only one visibility point at 110m. However, thanks to the 
spectral energy distribution of FU Ori, we concluded that 
the data were compatible with a standard accretion disk 
model. The binary scenario was also compatible with the 
data, but still required an accretion disk to explain the 
SED. 

We see here in the data the conjunction of 2 phenom¬ 
ena: a global decrease of the visibility with baseline and an 


The total complex visibility which results from these two 
objects is co mposed of the addition of the two terms (see 
lHodenll2nn(i for the mathematical details). 

5.1. Small-scale and bright structure 

In this section we focus our attention on the visibility 
behavior without considering the oscillations. Though we 
perform the model fits with all individual measurements, 
for clarity we plot the results with baseline-averaged visi¬ 
bilities, as displayed in Table El 


5.1.1. Fitting strategy 

We use a disk model and a methodology similar to those 
of paper I, but the data and technique used have under¬ 
gone several changes: the method to compute visibilities is 
now based on Hankel transforms, valid for axisymmetric 
distributions (See Eq. EJ ; in the disk model the exponent 
of the temperature profile power-law is a free parameter, 
allowing us to describe flared irradiated disks as well as 
self-heated viscous ones; visibility data taken with all base¬ 
lines are included in the study: the S EP now comprises 
observations of the lUezari et al.l ( 1 999l 5th edition 1 cata¬ 
log and from the 2MASS survey TCiitri et ^l2^^3^ . This 
model, although very simple, is able to reproduce a variety 
of situations without introducing new free parameters. 

The model consists of a geometrically flat and optically 
thick disk with an inner radius Tmin and outer radius r„iax, 
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100 

T—1—1—1—1—1—1—1—1—1—1—1—r 
X. = 1.63 microns 

PTI/NS 

T—1—1—1—1—1—1—1—1—1—1—1— 

X= 2.15 microns 

PTi/NS 


IOTA/S15N35 

^ VLTI/U1-U3 

IOTA/S15N35 

*> 

£ 0 

IOTA/S15N15 PTI/NW 

IOTA/S15N15 PTI/NW 

> 

# — 

< 

-100 

Ij_I_I_I_1_1_1_1_1_1_1_1_lI 

pirraw 

Ij _1_1_1_1_1_1_1_1_1_1_1_il 


-100 0 100 -100 0 100 


u(m) 


Fig. 1. (u,v) tracks corresponding to the observations of FU Orionis. Left panel: H data; right panel: K data. The 
(u, v) points have been labeled with the interferometer and the baseline used. 



Fig. 2. Calibrated squared visibilities of FU Orionis as function of hour angle in H (left column) and K band (right 
column). Rows from top to bottom display data obtained with the 6 interferometric baselines. 
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and features a radial temperature profile given by 


TW = T„(^) 


( 1 ) 


where tq = 1 AU is the reference radius, Tq the effective 
temperature at tq, and q a parameter^ usually ranging 
from 0.50 (flared irradiated disks) to 0.75 (standard vis¬ 
cous disks or flat irradiated disks). Each part of the sur¬ 
face of the disk located at radius r is assumed to emit 
as a blackbody at temperature r(r’), so that the flux and 
visibility are deduced by radial integration for a face-on 
disk with an inclination angle i = 0: 


FVO) = I 


rBx[T{r)] dr 


( 2 ) 


Vx{Bp,0) = 


1 2tt 

Fx{0)^, 


rBx[Tir)]Jo 


—R a 

X 


dr 


(3) 


with Jo the zeroth-order Bessel function of the first kind, 
Bx the Planck function, d the distance of the system and 
Bp a projected baseline. For a disk presenting an inclina¬ 
tion i and a position angle 0, these quantities become 


Fx{i)=Fx{d)cosR (4) 

EA(Bp,z)=EA(Bp(f,0),O), (5) 

where the corresponding baseline is 

S) = \l+ Bi g cos 2 i ( 6 ) 


and the coordinates are expressed in the reference frame 
rotated by the angle 9, 


Bu,g = Bu cos 9 — By sin0 (7) 

Byfi = By sin0 -I- By cos9. (8) 

Th e companion, FU Ori S, observed by IWang et al.l 
ll2nn4li and the contribution of the FU Ori central stellar 
source have negligible and compensating effects on the 
measured visibilities. FU Ori S is located 0.5arcsec away 
from the bright object, far outside the interferometric Held 
of view (viz. the “delay beam”, 9^ = ~ O.larcsec), 

and will therefore only combine incoherently with the flux 
from FU Ori N. The bias introduced by this incoherent 
flux corresponds to the contribution of this companion flux 
to the total, i.e. 1.6% decrease in visibility. In addition, 
we neglect the contribution of the central star, since it is 
unresolved and will only slightly increase the visibilities. 
The contribution is of the order of 2%, almost cancelling 
the decrease of visibility due to FU Ori S. 


5.1.2. Best fit and parameter constraints 

When trying to interpret the measured visibilities and the 
SED for the small-scale structure with the previous sim¬ 
ple disk model (neglecting the spot), the best fit has been 

^ See iMalbet fc BertouH ^199S^ ; iLachaume et alJ ll200ljl for 
details on the morphology of disks and their relationship with 
the q parameter. 


Table 4. Model parameters for the small-scale structure. 
This best fit has a reduced of 1.16 (with 306 individual 
measurements) and a “goodness of fit” of 2%. The 3-cr un¬ 
certainties on each parameter have been determined from 
a map around the minimum, as shown in left part of 
Fig. 0 


Accretion disc 


r ■ "i ® Rr^ r 

'min '-^•'^ — 1.8 'max 

To 7451^1 K q 

i 55^7 deg 9 


100 AU (fixed) 

f) 71+0.05 
C- ‘ J^-0.04 

47+( deg 


unambiguously found using a x^-gradient search method. 
The reduced x^ is 1.16, and the probability that statistical 
deviation from the model results in as large a deviation 
is 2% (“goodness of fit”). This implies that the model is 
unlikely to explain the observations (since we have been 
quite conservative in our determination of the data un¬ 
certainty), but we cannot rule out this possibility. This 
low goodness gives credence to trying the more compli¬ 
cated model with a spot (see Sect. 15.211 . The parameters 
obtained are presented in Tabled) with the standard un¬ 
certainty on the parameters. These errors have been de¬ 
termined from the x^ distribution displayed in left part 
of Fig. 0 — assuming a Gaussian error distribution. We 
checked that the uncertainties derived have values close 
to those predicted from the quadratic estimate of the x^ 
near its minimum. 

Unlike interfe rometric observations of T Tauri and 
Herb i g Ae stars llMillan-Gabet et ^ Il999l: lAkeson et alJ 
1200(1 ]Sio3 lEisner et al.ll2003[ l2004^ . the visibilities for 
FU Ori are compatible with the simple power-law radial 
temperature profile predicted by the standard disk model. 
We derive a temperature exponent q « 0.71 ± 0.05, which 
is actually strongly constrained by the slope of near-IR 
SED, and fairly consistent with the visibilities. The de¬ 
rived high temperature of about 745 ± 24 K at 1 AU hints 
to a viscous, self-heated disk rather than an irradiated 
disk. 

The shape of the visible SED sets a constraint on the 
cut-off radius rmin ~ 3-8 Rq. We find that the intensity 
of the near-IR emission strongly determines Tq cosi, while 
individual constraints on Tq and i are provided by the 
visibilities. The latter indeed determine Tq, rmin, * and 
9 (they have an influence resp. on the spatial extent of 
the disk in H and K, on the amount of central unresolved 
flux, and on the spatial extent in a particular direction), 
but their impacts are coupled. However, the visibilities 
also constrain the radial brightness distribution, provided 
that they be taken across a range o f wavelengths, as e x- 
plained by iMalbet fc Bergeil |2 to 2); iLachaun^ (|20^^. 
The visibilities at different baseline orientations constrain 

® These authors find a value for the temperature exponent 
q « 0.6 ±0.1 with a Icr uncertainty, which is different from the 
value derived in these papers. This discrepancy seems to come 
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X= 1.63 ixrni 

X= 2.15 |xm 

T 

lOTA/S 

15N15 


1 1 1 1 

lOTA/S 

15N35 



PT 

1 

CO 

. g. 



PTI 

NW 



PTI 

SW 



VLTI/L 

J1-U3 


ftoiv angle (JM) 


Fig. 3. Accretion disk model with baseline-averaged data corresponding to parameters listed in Table 01 Lines cor¬ 
respond to the model and symbols to measurements. The dotted line corresponds to an unresolved object. For the 
sake of clarity, data have been binned in the figure, but the model fit was performed using individual measurements. 
Top-left panel: spectral energy distribution and best fit model. Right panel: visibility data versus hour angle in H 
and K. Bottom-left panel: synthetic images in H and K in logarithmic scale. East is left and North is up. 
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Fig. 4. Disk-spot model corresponding to parameters of Table 0 Lines corresponds to the best model fit and markers 
to visibilities from Fig. [3 in binned form. Top-left panel: spectral energy distribution and best fit model. The solid 
line stands for the whole model and the dashed one for the spot. Right panel: visibility data vs. hour angle for each 
baseline in H and K. The contribution of the disk is displayed in dashed lines. Bottom-left panel: synthetic image 
in logarithmic scale. East is left and North is up. 
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cos(i) 9spot(°) 

Fig. 5. distribution of the disk alone (left column) and disk-spot (right column) models around the best fit (marked 
with a square). The white solid line draws the 3 -ct limit (“goodness-of-fit” > 0.25%). Top: inclination vs. temperature. 
Middle: temperature law exponent vs. minimum radius. Bottom left: inclination vs. position angle of the disk. 
Bottom right: position angle of the spot vs. projected distance to the centre. Left subpanels: x^ determined with 
both visibility and SED data. Middle subpanels: x^ determined with the SED only, right subpanel: x^ determined 
with the visibilities only. 


the inclination and the position angle of the disk. The 
visibilities at short baselines put no constraints on the ge- 
ometry of the dis k because it is only marginally resolved 
ijTjachaiimelbnnfiji . but we can rule out the contribution 
of an extended (sub-arcsec) structure (e.g. scattering by 
an envelope) because IOTA visibilities at < 30 m are 
unresolved'*’. 

Figure 13 illustrates the comparison between the model 
and the observations. The SED data and the result of the 
fit are displayed in the upper left panel. The visibilities 
averaged by baseline, listed in Table |3 are presented in 
the right panel with the visibility curves corresponding to 
the (u, v) track of each baseline. Finally in the lower left 
panel, we show the corresponding synthetic image of the 
disk. 

5.2. Large-scale structure 

In this section we focus our attention on the high- 
frequency oscillation in the visibility function, mostly vis¬ 
ible in the 110m AT-band data of PTI. We believe it to 
come from an off-center unresolved light source that we 

from the reduced {u,v) coverage used in that paper and from 
their ignoring the SED (which pushes q towards 0.75). 

'* An extended contribution yields a quick decrease of the 
visibilities with baseline. 


call a “spot”. We fit to the data a compound source model 
consisting of an accretion disk (as described above) and 
an unresolved spot. The latter is modeled as a circular 
black-body of uniform brightness with a temperature Tspot 
and a radius rspot- The two components are separated by 
the projected physical distance dspot with a position angle 
fi*spot (with a 180deg ambiguity). The fitting procedure is 
the same as in Sect. o 

5.2.1. Properties of the spot 

The best model fit was obtained with the parameters listed 
in Table |3 The reduced x^ is 0.95 and the probability 
that statistical deviation from the model accounts for as 
large a deviation is 70%, so the fit appears as a consistent 
explanation for the data - in comparison to the 2% of the 
disk-alone model - but the detection of the structure is 
only at a 2-tT level. The best-model fit is shown in Fig. 01 
and the probability distribution around the minimum is 
shown in the right part of Fig. [3 

The disk parameters and their uncertainties are con¬ 
sistent with the results of Sect. 15.ll except for the position 
angle: its uncertainty is much larger in our last model fit 
(8 ± 20deg) than in the previous one (47 ± 10 deg), and 
both values differ by 2a. This angle is constrained (to¬ 
gether with the inclination) by the visibility difference at 
three different baseline angles, namely with the NW, SW, 
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Table 5. Model parameters for the large-scale structure 
without secular variation. The reduced of the best fit is 
0.95 (with 306 individual measurements) and the “good¬ 
ness of fit” is 70%. 



Binary system 


I a 

Q-spot 

AU 

n a 

"spot 

-50t)deg 


Accretion disc 


rmin 

4.91^;? Re 

^^max 

100 AU (fixed) 

To 

692+2® ^ 

q 

o.7i^°:®® 

i 

48 tgdeg 

9^ 

si®) deg 


Uniform 

stellar disc 

rji b 

spot 

4600l»°° K 

"^spot 

5.0 Re (fixed) 


Other distant minima exist. 


^ Quadratic error from the fitting routine 


and NS PTI baselines (see the right panel of Fig.0J. When 
a spot is added, this difference is changed, resulting in a 
different position angle: in our best fit, the disk visibility 
for the NS baseline is unchanged since it is the average 
of the oscillations, while it is higher for the NW and SW 
baselines because the measured visibility (disk plus spot) 
is in the low of the oscillations. 

As Fig. 0 shows (see bottom panel, Eis a function 
of the spot location), there are multiple acceptable “best” 
model fits featuring radically different spot locations. The 
aforementioned solution is actually the best in terms of 
modelling the PTI/NW and NS oscillation in K, but was 
not found to be statistically better than others that do 
not feature such oscillation. We think, however, that the 
oscillation is actual and that the model ambiguity results 
from the poor accuracy of other measurements (IOTA in 
H and K, and PTI in H). Thus, results on the spot location 
should be taken with caution. 

We have chosen to parametrize the flux of this spot 
in terms of stellar photosphere with the temperature and 
the radius of the photosphere as parameters. However only 
the K data can be used to constrain them and we have 
fixed the value of the radius to 5Rq, which is typical for 
a young star. The derived temperature depends strongly 
on the value taken for the radius and therefore should be 
considered with caution. 

The location of the spot is quite well constrained by 
the shape of the oscillation with a projected distance to 
the center of 10 ± 1AU and a position angle —50 ± 1 deg 
(with 180 deg ambiguity). The amplitude of the oscilla¬ 
tions in the K-band constrains the flux ratio of the spot 
to the primary, AK « 3.9 ± 0.2 and Aff « 3.6 ± 0.2 mag. 
However, there are other acceptable minima (goodness- 
of-fit > 50%, see the bottom right x^ maps in Fig. 
which do not reproduce the oscillation, but they cannot 
be ruled out. We should note that the oscillations are de¬ 
tected mostly in the K band and therefore AJI is not 


Table 6. Parameters of the best model fit of the location 
of the spot in 1998, 1999, 2000, and 2003, and its quality. 


year 

4pot (AU) 

6'spot (deg) 

X® (total) 


1998 

8.0 ±1.9 

-48.0 ±3.3 

0.63 

0.54 

1999 

9.6 ±1.5 

-49.4 ±2.2 

0.95 

0.78 

2000 

10.1 ±0.6 

-50.3 ± 1.6 

1.04 

0.91 

2003 

11.1 ±0.3 

-49.9 ± 1.7 

2.54 

4.02 


constrained by both the oscillations and the model, since 
there is little chance that a spot at 4000 — 5000 K be a 
pure black-body. 

5.2.2. Secular evolution 

Assuming that the spot location of the disk-spot model 
is that derived in Sect. 15.2.11 (which allows us to explain 
the PTI oscillation in K), we carried out a model fit of 
its location for each year of observation. The results are 
displayed in Table 0 The fit is successful in the years 
1998, 1999, and 2000 and feature a 3-a detection of an 
almost radial motion. The data obtained in 2003 are not 
sufficient to constrain the position of the spot for this year 
explaining the bad 

6. Discussion 

6.1. Accretion disk versus photosphere 

The current data set confirms the result from paper I and 
shows consistent visibility data. We can definitively con¬ 
firm that we are observing an object whose size (Gaussian 
full width at half maximum) is 1.5mas in K wide corre¬ 
sponding to 0.7AU or 150i?@ a t a distance of 450 pc (see 
paper I for details). Therefore if lHerbig et all ll2003li is cor¬ 
rect in interpreting the high resolution spectroscopic data 
as stellar chromospheric activity, then the star must be ac¬ 
companied by cooler circumstellar material that spreads 
beyond 15Oi?0 that is probably a circumstellar disk. 

The temperature Tq obtained in our fits yield an ac¬ 
cretion rate® of M = (6.3±0.6) x lO“®(M,t/M 0 )“^ MQ/yr 
for the disk alone and of M = (4.7 ± 0.6) x 

1 O“®(M*/M 0 )“^ M 0 /yr for the disk accompanied by 
a spot. The effective temperature-accretion rate rela- 
tion is the one of the st andard model (see Eq. 2.7 in 
IShakiira fc Sunvae^ll973^ . It is interesting to note that 
the corresponding temperature at the inner radius is 
10050 K (resp. 10140 K). In our present model, the disk is 
optically thick whatever the physical process which dom¬ 
inates the opacity (gas or dust). A more detailed study of 
this boundary region would be of great interest. 

® Uncertainties include the fitting uncertainty and the uncer¬ 
tainty on the boundary condition in the standard model (given 
by the stellar radius, liberally assumed to be between 0.5 to 6 
solar radii). 
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Compared to paper I, we are able to begin constraining 
the geometry of the accretion disk, with an inclination 
angle of the order of 50 degrees and position angle of the 
order 10-40 degrees. More accurate data and better [u, v) 
coverage is still desirable, but we already show the path 
toward such constraints. 

Other observa tio ns of Herbig Ae stars 
( Millan-Gabet, et ^ I 2 OOII: lAkeson et al.1 l2002t 

lEisner et all l2003l l2004ll and o f T Tauri stars 
I Akeson et al.l l20^ IColavita et al.l bOOJljl cannot be 
interpreted in the framework of the standard disk 
model. Measured near-IR visibilities are smaller than 
the theoretically expected values, suggesting puffed-up 
inner disk walls, and flaring is often necessary to fit 
mid- and far-IR ph otometry jPullemond et all l200lt 
iMuzerolle et al.ir200.lll . In contrast, our wor k on FU Ori, 
as wel l as some Herbig Be stars observed by lEisner et al.l 
1 I 2 OO 4 I 1 . shows that the observed circumstellar disk 
structure is consistent with the standard accretion disk 
model, i.e. geometrically flat and optically thick with a 
temperature law index close to 0.75. We suggest that 
accretion-dominated objects (such as FUors and Herbig 
Be stars) fit the standard accretion disk model, while 
Herbig Ae stars and T Tauris may have different disk 
structures due to the influence of irradiation. This is also 
consistent with the fact that in FUors the disk completely 
dominates the luminosity and therefore the physics of the 
disk would be expected to be governed by accretion and 
no t by external stellar heating. 

iHartmann fc KenvonI (jl 985|i have presented an inter¬ 
esting interpretation of double lines o bserved in FU Ori by 
spectroscopy. The controversy with iHerbig et alJ lj200,‘l(l 
comes from the fact that these authors are not able 
to localize the origin of the se lines. The explanation of 
IHartmann fc KenvonI 1 I 1 985^ appears to be more consis¬ 
tent with our data. 

The interferome tric instrument AMBER on the VLTI 
llPetrov et al ]|20^ has the potential to detect these dou- 
ble lines and identi fy the region where they are emitted 
llMalbet et alJlioO^ . 

6.2. The nature of the bright spot 

Our analysis of the data shows that the scenario which 
involves a bright spot is more likely to explain our obser¬ 
vations. The presence of such a spot would modulate the 
visibility curves. Even if the amplitude of the oscillations 
is not large and the number of cycles is small, we believe 
that this spot is indeed present. As can be seen in Fig. 0] 
the ripples of the visibility also explain some low visibili¬ 
ties obtained at IOTA with short baselines. Thi s detection 
was already pointed out bv iBerger et alJ ll200fil . based on 
a smaller data set, and the latest observations confirm this 
behavior. In addition, our analysis by year, using indepen¬ 
dent sets of observations, yields parameters which are close 
to each other. To confirm the presence of this bright spot, 
we suggest observing FU Orionis in several wavelengths 


and possibly with spectroscopy (see Fig.^lfor an example 
of the correlated spectrum that would be obtained with 
a specific configuration of the VLTI with the instruments 
AMBER and MIDI). 

The nature of this bright spot is still unknown. This 
structure, assuming it is unresolved, could be due to an 
increa se of the intensity in the disk due to a thermal o ut- 
burst dClarke et al.lll99(l iBell fc LinlllQQ^ iBelllllQ^ . If 
the assumption of a photosphere radius of 5Rq is correct 
(see remark in Sect. KTU . the derived color temperature 
of 4600 K of this spot is similar to temperatures observed 
in stellar photospheres and therefore is consistent with the 
presence of a young stellar companion. 

The year-by-year analysis indicates that the spot may 
be moving with a velocity of about 1.2 ± 0.6 AU • yr“^ 
on a rectilinear trajectory. We have insufficient data 
to complete the analysis but additional data will pro¬ 
vide more time-sampling to confirm or exclude the mo¬ 
tion of the spot. The proper motion of FU Ori is 
116 ± 104 X 10~®s/yr in right ascension and 86.4 ± 
15.5 mas/yr in declination (Carlsberg Merid ia n Cat alogue: 
ICopenhagen University Observatory et~ ] lIlOOflM . i.e. 
39 ± 7AU/yr at position angle of 12 ± 10 deg. The spot 
motion, if actual, is therefore not due to a foreground or 
background star. It could be explained by the motion of 
a companion star orbiting on a very highly eccentric orbit 
around the primary star or an orbit inclined with respect 
to the disk. 

This explanation is consi s tent w ith the scenario intro- 
du ced bviBonnell fc Bpt ienI lll992ll and recently revisited 
bv lHeipurth fc Asninl l|2004bll following the dis covery of a 
binary companion t o FU Ori ||W anfi_etaI ■I 2 OO 4 II . The com¬ 
panion observed bv lWang et al.1 ll2004ll cann^ be respon¬ 
sible for or be the consequ e nce oL the 1936 outburst as 
noted bv lReipiirth fc AspinI I 2004bll . However these same 
authors have also pointed out that “FU Orionis itself must 
be a close binary, with a semi-major axis of 10 AU or less”. 
Confirmation and extension by future measurements are 
needed to confirm that the spot we are observing might 
indeed be the companion they were looking for located on 
an inclined or highly eccentric orbit. 

This spot could also represent t he signature of the pres- 
ence of a young planet. Recently llLodato fc Clarkell2004ll 
proposed that FUor outbursts could be caused by a planet 
at about lOi?© which causes gap instabilities. In the case 
of FU Ori, the mass of the planet would have had to be 
~ 15Mjupiter in order to explain the rapid rise of the lu¬ 
minosity up to 500Lq. We do not have enough color mea¬ 
surements to derive the temperature and the luminosity 
of the spot in order to check that it is compatible with 
a 15Mjupiter planet. Nor can we explain why the planet 
would be migrating away from the center. 


7. Conclusion 

We have presented the most complete set of interferomet¬ 
ric observations obtained so far in the domain of star for¬ 
mation with several infrared interferometers: PTI, IOTA 
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Fig. 6. Predicted squared visibility curve of the disk+spot 
model for FU Ori through the infrared spectrum with the 
triplet of VLTI 8 m telescopes U2-U3-U4. 

and VLTI. These observations help us to better under¬ 
stand the exact nature of FU Orionis. Our measurements 
confirm that FU Ori hosts an active accretion disk with a 
« 5 X lO“®M 0 /yr accretion rate. We have also marginally 
detected an unresolved bright structure, that we identify 
as a possible close companion on a highly eccentric orbit. 
This putative companion may be responsible for the out¬ 
burst of the accretion rate that led to the rapid increase 
of the luminosity. 

These observations also show that the inner structure 
of YSOs is rather complex and cannot be fully understood 
with only a small number of interferometric observations. 
More wavelength and (m, v ) coverage is required to fully 
constrain models. The next step will probably be the re¬ 
construction of an image of this interesting object based 
on extensive observations with the VLTI. 
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